Shape coexistence and tilted-axis rotation in neutron-rich hafnium isotopes 
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We have performed tilted-axis-cranked Hartree-Fock-Bogoliubov calculations for a neutron-rich 
hafnium isotope (^*^Hf) whose proton and neutron numbers are both in the upper shell region. We 
study whether the shell effects play a role in producing high-_K' isomers or highly gamma-deformed 
states at high spin. In particular, the possibility of shape coexistence and the effect of wobbling 
motion are discussed. 
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Hafnium isotopes {Z — 72) are best known as nuclei 
that have high-X isomers (e.g., the = 16+ isomer 
in i^^Hf, with a half life ti/2 = 31 yr §). From a 
viewpoint of the Nilsson model, a reason is that proton 
single-particle levels are filled up to the upper part of 
the shell where there are many high-fi states at prolate 
deformation (f2 is the angular momentum projection on 
the nuclear symmetry axis). The presence of long-lived 
high- if isomers indicates the existence of axial symme- 
try to make K = J2i a good quantum number. With 
recent developments of experimental techniques, such as 
fragmentation |Q and deep-inelastic reactions |^ in pop- 
ulating high-spin states, the study of nuclei in the well 
deformed rare-earth region moves away from the P sta- 
ble line towards the neutron-rich isotopes. D'Alarcao 
et al. recently discovered several high- if isomers in the 
neutron-rich hafnium isotope of ^^^Hf Q. For neutron- 
rich hafnium isotopes, whose neutron Fermi level is lo- 
cated in a similar position to the proton one (i.e., the 
upper half of the shell), we can expect an even more im- 
portant role for high- if isomeric states At the same 
time, however, the empty nucleon states near the Fermi 
surfaces can be considered as hole states, and these may 
induce substantial gamma deformation which breaks the 
axial symmetry. 

A few theoretical investigations of high angular mo- 
mentum collective states have been reported for such 
neutron-rich hafnium isotopes: by means of a microscopic 
method (cranked Hartree-Fock-Bogoliubov method, or 
C-HFB), the possibility of a collective oblate deforma- 
tion at high-spin (i = 26?i) in ^^"Hf was predicted 
by means of a macroscopic-microscopic method (total 
routhian surface calculations, or TRS), the existence of 
non-collective prolate-deformed states becoming yrast, 
i.e., lowest in energy, for I > lOTi was shown in i82-i86jjj: 
0. In the latter study, the collective oblate deformed 
states become yrast when quite high in spin (i = 36h in 
^^^Hf). In both of the studies, the rotational states are 
treated only by the one-dimensional cranking model, but 
recent developments in the tilted-axis-cranking model 
|p[-p^ allow for the possibility of further investigation 



of these neutron-rich hafnium isotopes. 

In this paper, in a microscopic framework of the two- 
dimensional tilted-axis-cranked HFB method, wc analyze 
the high-spin structures predicted for the neutron-rich 
hafnium isotope ^^^Hf (which has N/Z = 1.53). 

First, the calculational procedure is briefly described. 
We solve the HFB equations with a 2d-cranking term self- 
consistently for ^*^Hf, following the method of steepest- 
descent Nucleon numbers and the total angular mo- 
mentum are constrained during the iterations; (Nt) — 
Nr (r = p,n); = (4) = J cos 6*; Jy = (4) = 0; J, = 
(iz) = Jsm9. (Note that the tilt angle 9 is measured 
from the a;-axis to see the deviation of the rotation axis 
from the x-axis. In this study, it corresponds to a princi- 
pal axis of the quadrupole deformation, perpendicular to 
the symmetry axis that is chosen to be the z-axis). We 
also constrain off-diagonal components of the quadrupole 
tensor such that the intrinsic coordinate axes coincide 
with the principal axes of the quadrupole deformation. 
See Ref. for details. 

The one-dimensional cranking calculations (or Id- 
cranking calculations), in which the rotation axis is fixed 
to be along the x-axis (6 = 0°), are performed as follows. 
First for J = 0, the Nilsson -I- BCS state is taken as a 
trial state. Then cranked HFB states, ^'u(J), are calcu- 
lated up to J = 40?i with increments AJ = O.lTi. This 
way of calculation is what we call "up-cranking" calcula- 
tions. Then, using the up-cranking solution at J = AOh 
as a trial state, "down-cranking" calculations are simi- 
larly performed from J — iOh down to J = to obtain 
*d(J). 

The calculations of tilted-axis-cranked HFB states (2d- 
cranked states, or ^*( J, 9)) are carried out by starting at 
a state ^''( J, 9 = 0°), which are calculated through an up- 
cranking calculation, and performing a "forward" tilting 
calculation up to = 90° with increments A6' = 0.5° 
for each (integer) J. Then in a similar manner we make 
"backward" tilting calculations from 9 = 90° to = 0°. 
At certain tilt angles, the forward and backward tilting 
results do show interesting differences, particularly, at 
high spins. 
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Our Hamiltonian consists of a spherical part (Hq) and 
a residual part (Kos) employing the pairing-plus-Q • Q 
interaction. In the hafnium isotopes the hexadecapole 
interaction can be important, but we have checked that 
our results below are not affected very much by including 
the interaction. Thus we omit it in the present study. 
The Hamiltonian is thus written, 
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in which means a spherical Nilsson level and i runs 
all over the model space. Force parameters n and Gr 
are determined in the framework of the Nilsson + BCS 
model by giving input parameters for the quadrupole de- 
formation (/?'"', 7™') and gap energy (A™'). (Our defini- 
tion of 7 is taken from pp.8 in Ref. and is opposite 
to the Lund convention in the sign.) In this paper, we 
employ a set of the input parameters based on the cal- 
culations by MoUer et al. |Q because we find it gives a 
good agreement with experimental data. For ^^^Hf, we 
employ A|f' = 0.725 MeV; Ajj" = 0.625 MeV; 7''^' = 0°; 

= 0.270 or 0.268. We will see soon the reason for the 
two values of 

Our single-particle model space is almost the same as 
the choice of Kumar and Baranger (two major shells 
in the spherical Nilsson model: = 4, 5 for proton and 
iV = 5, 6 for neutron), with two extra single-particle or- 
bits, proton ii3/2 and neutron ji5/2- The single particle 
energies are the spherical Nilsson model energies with 
A-dependent Nilsson parameters jlj] . 

First of all, let us discuss a principal axis rotation 
(6 = 0°) in Id- and 2d-cranking calculations. Fig.l shows 
the energy spectra of one-dimensional up- and down- 
cranking states ( 'I'u(>/) and 5'd( J), respectively ) and the 
states oi 9 — 0° in the 2d-cranked calculations,^* (J, 9 = 
0°),that are obtained through the backward tilting pro- 
cedure at a given value of J. They are plotted also in 
the inset of Fig. 1 together with the known experimental 
values of the g-band. The g-band is well reproduced by 
our one-dimensional cranking calculations. 

We have two kinds of Id-cranked states plotted in 
Fig.l. They are calculated with the same conditions 
except for the initial value of /3'"'. One (denoted as 
"A" in the graph) is obtained with the initial value of 
/3™' = 0.270 while the other (denoted as "B") is with 
/3'"' = 0.268. The small difference in gives almost no 
difference in the wave functions at low spin (J < IQh), 
but it does lead to a significant difference at high spins. 
The sensitive dependence of the high-spin HFB solutions 
to shows that the energy manifold in the variational 
space near the crossing regions has several local minima 
that are almost degenerate. The states "A" and "B" im- 



ply the possibility of band crossings 0, which correspond 
to the regions Uh < J < 26h and at 22h < J < 2,lh 
respectively [| in the calculations, but a substantial dif- 
ference between the states is seen in the gamma deforma- 
tion. Before the crossing regions both of the states have 
7 ~ 10°, while after the crossing regions the states "A" 
have near-prolate deformation with a negative gamma 
value (7 ~ —10°) and the states "B" have oblate shape 
(7 ~ 60°). (See Fig. 2(a) for evolutions in gamma de- 
formation for each solution). From an analysis of our 
numerical results, the oblate deformation is caused by 
the gradual alignments of both ii3/2 neutrons and hii/2 
protons, while, in addition to these alignments, the near- 
prolate deformation with 7 ~ —10° is caused by the quick 
alignment of j25/2 neutrons. (Note that the neutrons in 
the ji5/2 orbits are not part of the usual Kumar-Baranger 
model space). The solution "A" is reported in this pa- 
per for the first time. Xu et al. found that there is no 
stable minimum corresponding to this solution in their 
TRS calculations 0|. However, they performed their 
(Id) cranking calculations for given rotational frequen- 
cies {ijJx), while we have performed calculations for given 
(average) angular momentum vectors (J^, Jy, Jz)- 

In a spin region 23?l < J < 35?i, the three kinds of 
states having different gamma deformations (7 ~ ±10° 
and 60°) are close to each other in energy. This result im- 
plies a manifestation of shape coexistence at high spin, 
or multi-band crossings among bands specified by dif- 
ferent gamma deformations (or corresponding rotational 
alignments) . 

It is interesting to see the energies for the state 
\I'*(J, 6* = 0°), which are represented by the thick solid 
line in Fig.l. The line has two discontinuities, at J = 27?i 
and 34?i, implying two configuration changes. At lower 
spin (J < 26?i), «-*( J, 6* = 0°) follows the Id-cranking cal- 
culations which give rise to nearly prolate shape (7 ~ 10° 
at J ~ 25?i). Then at J = 27?i, the gamma deformation 
changes to 7 ~ —10°, which is the same as "A". Finally, 
at J = 34S,, 5'*(J, 6* = 0°) changes into the same state as 
"B", having oblate shape (7 ~ 60°). The result that the 
states having different gamma deformation are connected 
by the tilted-axis-cranking solutions indicates the impor- 



*As Hamamoto et al. pointed out the validity of the 
cranking model is questionable in the band crossing region be- 
cause of the semi-classical aspects in the model. However, the 
model works well outside the crossing region and could explain 
qualitative features near the crossing region, such as align- 
ments. For a more acculate analysis of the crossing region, 
we should employ the so-called "variation-after-projection" 
method jll], the generator coordinate method jisj, or the 
diabatic method |l6j. 

^Note that we do not mean here that these values of the 
crossing angular momentum correspond exactly to the exper- 
imental values. In general, the simple self-consistent cranking 
calculations do not reproduce the value precisely 11 7|. 
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tance of the tilting degree of freedom [9) for the search 
for excited states near band crossings. In Fig. 2(b), we 
show how the quadrupole deformations (/3 and 7) evolve 
at J = 34?i as we vary the tilt angle {6). The corre- 
sponding energy curves are plotted in Fig. 3 (a). There 
are three types of solutions with different gamma defor- 
mation: (i)the solution having no ^-dependence in tilt 
angles for < 20° is oblate (7 ~ 60°); (ii)the solution 
which shows a tilted rotation minimum at ~ 10° has 
negative gamma deformation (7 ~ —10°); (iii)the solu- 
tion with a minimum at 9 = 90°, which may correspond 
to high-K states, has 0° < 7 < 10°. 

The energy difference between the states of type (i) and 
(ii) is roughly constant and small (~ 500 keV), so that 
these two states can couple to form states with mixed 
deformation. The energy curves for these states are shal- 
low in the range 0° < 9 < 30°, so that fluctuations in the 
rotation axis, or wobbling motion jist , can be expected. 
However, gamma deformations for each state of type (i) 
and (ii) are quite constant against variation in the tilt 
angle up ~ 30°. Therefore, rather than a picture in 
which states of type (i) and (ii) are coupled thorough the 
wobbling motion, we should have a picture where they 
are connected possibly through 7 tunnelling, and where 
the mixed states wobble around the tilted rotation mini- 
mum at ~ 10°. Nevertheless, a coupling of these mixed 
states with states of type (iii), is possible through wob- 
bling motion. 

It is interesting to see in Fig. 3(a) the energy curve at 
J — AOTi corresponding to the states of type (i) above. 
There are two minima at = 0° and 90°, but the barrier 
height between them is only ~ 50 keV, not visible in the 
plot. The corresponding deformation is collective oblate 
(the y-axis is the symmetry axis) and almost constant 
over the entire range of 9. For a strict oblate symmetry, 
there would be no energy dependence on 9. The projec- 
tion of the nuclear deformation onto the x — z plane is 
a circle, so that there is no preference for a direction of 
(collective) rotation in the x — z plane. For finite triaxi- 
ality, it is important to consider wobbling motion in the 
9 direction. 

From these discussions, we can deduce that the 
shape coexistence creates successive backbends (sudden 
changes in moments of inertia) in the excited rotational 
bands of ~ 0°. The first backbend, which is caused by 
neutron j]^5/2 alignment, in addition to alignments of neu- 
tron ii3/2 and proton hii/2: is expected at lower spins, 
which corresponds to 17?i ^ J ^ 26h in our calcula- 
tions. The second one, which is caused by de-alignment 
of the ji5/2 neutrons and retaining ahgnments of the ii3/2 
neutrons and hj^i/2 protons, can be seen at higher spins 
(corresponding to 22?i ^ J ^ 31?i in the calculations). 
According to Fig.l, the second backbend can be more 
pronounced than the first one. It is also possible that 
these two backbends are mixed together to create one 
giant backbend as Hilton and Mang predicted but 
according to our analysis three types of gamma deforma- 



tions are involved: 7 ~ ±10° and ~ 60°. 

Now, let us look at the calculated high-X states. In 
Fig.l, states of 6I = 90°, **(J,6l = 90°), are shown with 
symbols. These states correspond to a local mini- 
mum a.t 9 — 90° in the energy curve (see Fig. 3(b)). In our 
calculations, this minimum starts to appear at J ~ 8h, 
and becomes the lowest minimum at J — 12h and at 
higher spins. The 9 — 90° minimum may be considered 
approximately to correspond to a high-if state, and the 
corresponding 2-axis cranked state, or ^'*( J, 9 = 90°), is a 
simulation of nuclear rotation where single-particle angu- 
lar momenta carry most of the total angular momentum. 
We have checked that the corresponding gamma defor- 
mation is almost zero as in Fig. 2(b) which shows that 
there is axial symmetry for these states. 

However, we should note that Jz{= J sin 6) in the self- 
consistent cranking calculations is not an eigenvalue but 
just an expectation value, due to the rotational symmetry 
breaking by the mean field. Besides, in the tilted-axis- 
cranked HFB calculations, angular momenta are fully 
mixed in the sense that even a mixture among even and 
odd angular momenta happens. This is because the sig- 
nature symmetry, a discrete subgroup (-D2) of the rota- 
tional group 0(3), is broken by tilted rotation |T^. We 
should therefore keep in mind that the mean-field de- 
scription of high-X states has a certain limitation. 

Remembering the above remarks, let us look at our 
results for 5"* (J, 6* = 90°). The experimental energy of 
the isomer, which is tentatively assigned to = 13+, is 
2.572 MeV relative to the ground state energy, while the 
numerical values are 2.101 MeV (J = 12?i), 2.255 MeV 
(J = 13n), and 2.408 MeV (J = 14?i). Deformations of 
*'(J = 13h,9 = 90°), are calculated to be /3 = 0.2661 
and 7 = 0.015°. 

We can consider whether this isomer is yrast or not. 
(^Experimentally, the g-band is identified only up to 8?i 
[|l[). Xu et al. calculated that the non-collective prolate 
state would be yrast at J = 13?i 0. Our self-consistent 
calculations also show that 5'*(J, 6* = 90°) is lower in 
energy than principal axis rotation, '^*{J,9 — 0°), for 
12?i < J < 34?i, so that in this region the high- A' isomeric 
states can be favoured relative to the collective rotation. 

Let us use the term "tilted rotation" for the states with 
local minima, with 6* ^ 0, 90°. At J = 13?i (see Fig.3(b)), 
there is tilted rotation {9 ~ 15°) with an excitation en- 
ergy of 2.593 MeV, which happens to be quite close to 
the experimental value of the K'^ = 13+ isomer. Numer- 
ical deformation values for this minimum are /3 — 0.2704 
and 7 — 4.62°, that is, the shape is nearly prolate. It is 
possible to consider the tilted rotation to describe a rota- 
tional member of a high- A' state. However, at this angu- 
lar momentum, the potential energy curve is shallow in 
the whole range of 9 (see Fig. 3(b)), and the barrier height 
is the same order of magnitude as uncertainties from the 
mean-field approximation. The question as to whether 
the high- AT state should be described as either a tilted 
rotation state or a z-axis cranked state 5'*(J, 61 — 90°) 
(or coupling of them) should be answered by a quantum 
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mechanical calculation by using angular momentum pro- 
jection (and/or the generator coordinate method), which 
we plan to study in the future. At this moment, the 
description of high-ii' states in the framework of self- 
consistent tilted-axis-cranking calculations is reasonably 
good up to an accuracy of several hundred keV, but more 
detailed studies are surely necessary. 

For these shallow energy curves, we can consider the 
wobbling motion to relate and cause transitions between 
high-ZsT isomers and low-K states such as the g-band. For 
the energy curve at J = 13h, there are three minima at 
= 0° (principal axis rotation), 9 = 15° (tilted rotation), 
and e = 90° (possible high- if states) within 600 keV. If 
the first minimum represents a rotational member in the 
g-band and the third (and/or second) minimum repre- 
sents the one in a high-if band, then the correspond- 
ing band crossing is expected to show (experimentally) a 
strong coupling as evidence for a realization of the wob- 
bling motion. 

In summary, we have performed tilted-axis-cranked 
HFB calculations for ^®^Hf, and investigated high-spin 
states near the yrast line up to J = 40h. For our param- 
eter set, the comparison of the experimental data for the 
g-band and a high- if isomeric state {K^ = 13+) with 
our calculation gives a reasonable agreement within the 
framework of the mean-field approximation. With our 
modified single-particle model space based on Kumar- 
Baranger's choice, we found a new HFB solution with 
near-prolate deformation (7 ~ —10°) involving j 15/2 neu- 
tron alignment at high spin (J > 17fi). We discussed 
the possibility of shape coexistence among three states 
in principal axis rotation with different gamma deforma- 
tion: two with near-prolate shapes (7 ~ ±10°) and the 
other with oblate shape (7 ~ 60°). An analysis of the 
possible backbends in yrare states was also given. In ad- 
dition, we discussed the effect of wobbling motion as a 
coupling mode between low-i^" and high-7f states. 
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FIG. 1. Energy spectrum near the yrast line obtained through the tilted-axis-cranking calculations. In the large graph, we 
show the energy with an arbitrary subtraction of rotational energy, 0.005J(J + 1) McV for 10ft < J < 40?;-, in order to see 
excited structures in detail. In the small graph, the energy curve (without subtraction) is shown for the full range of angular 
momentum Oh < J < 40h. Experimental values for the g-band (open circle) and a = (IS"*") isomer (asterisk) are plotted for 
a comparison. Ono-dimonsional up- and down-cranking arc depicted by upward and downward triangles, respectively. When 
the up- and down-cranking results are the same, the triangles are superposed to give the stars. Two kinds of Id-cranked 
calculations are shown as "A" and "B", which are obtained by = 0.270 (open triangles) and 0.268 (closed triangles), 
respectively. "A-l-B" means that the two Id-cranking solutions "A" and "B" give similar results. Also "u-|-d" means that 
the up- and down-cranking calculations show no differences. The thick solid line indicates states of 6 = 0° obtained in the 
2d-up-cranking calculations (*!'*( J, ^ = 0°)) while plus signs denote 2d-up-cranking states ^^{J,6 = 90°). The tilted rotation 
minima are represented by diamonds (see Fig.3(b)). 
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FIG. 2. (a) Gamma deformation in Id- and 2d-cranking calculations are shown with respect to J(= Jx)- The symbols are 
the same as those in Fig.l. (b) Gamma and beta values at J = 34?i in the 2d-cranking calculations are shown with respect to 
6. "Fwd" means the forward tilting calculations, while "Bwd" means the backward calculations. 
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FIG. 3. Energy curves with respect to tilt angle 6. (a)energy curves at J = 34h and 40h. Three states having diflFerent 
gamma deformation are energetically close in the range 0° < 9 < 40° . For J = 40% there are two minima, at ^ = 0° and 90° , 
the barrier height between them being only about 50 keV. (b)energy curve at J = 13h. For 12h < J < 20h, there are three 
minima sd. 6 = 0° , 90°, and ~ 15°. The last minimum implies tilted rotation, (c)energy curve at J = 4h. In the low-spin region 
(J < 12?i), there is only one minimum, at 9 = 0°, i.e., principal axis rotation. 
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